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ABSTRACT 



We report the discovery of diffuse hard (1-5 keV) X-ray emission around 
Jupiter in a deep 160 ks Suzaku XIS data. The emission is distributed over 
~ 16 x 8 Jovian radius and spatially associated with the radiation belts and the Io 
Plasma Torus. It shows a flat power-law spectrum with a photon index of 1.4±0.2 
with the 1-5 keV X-ray luminosity of (3.3±0.5) xlO 15 erg s _1 . We discussed its 
origin and concluded that it seems to be truly diffuse, although a possibility of 
multiple background point sources can not be completely rejected with a limited 
angular resolution. If it is diffuse, the flat continuum indicates that X-rays arise 
by the non-thermal electrons in the radiation belts and/or the Io Plasma Torus. 
The synchrotron and bremsstrahlung models can be rejected from the necessary 
electron energy and X-ray spectral shape, respectively. The inverse-Compton 
scattering off solar photons by ultra-relativistic (several tens MeV) electrons can 
explain the energy and the spectrum but the necessary electron density is >10 
times larger than the value estimated from the empirical model of Jovian charge 
particles. 

Subject headings: planets and satellites: individual (Jupiter, Io) — X-rays: general 
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Introduction 



In situ measurements and ground-based observations have shown the existence of 
Jupiter's high-energy elect ron Van Allen radiation belts with evidence for electrons at 



energies of keV to 5 MeV (IBerge Gulkis 1976 



Bolton et al. 20021 ). Theoretica 



Van Allen et al. 1975 



Fischer et al. 1996 



via the inward radial diffusion ( Goertz et al. 1979 



ly, electrons are energized by the beta tron acceleration 



Miyoshi et al. 19991 ) and recently non 



adiabatic acceleration via wave-particle interactions has been discussed (IHorne et al. 20081 ). 



The trapped electrons lose their energy through s ynchrotron radio em ission and some 



fraction are absorbed by Jovian moons and rings (IBagenal et al. 20041 ). However, there 
remains much debate over the precise energy spectrum and spatial distribution of the 
energetic electrons in the radiation belts. 

Jupiter is the most luminous planet in the solar system at X-ray wavelength 



f lMetzger et al. 1983 



Waite et al. 1997 



Bhardwaj et al. 20071 ). It emits X-rays from mag- 
netic poles by the charge exchange interaction with the solar wind and magnetospheric heavy 



ions and also by the bremsstrahlung emission via energetic electrons ( [Gladstone et al. 2002 
Branduardi-Raymont et al. 2004 ; iBranduardi-Raymont et al. 2007af ). A typical X-ray 
luminosity of the auroral emission is ~10 16 erg s _1 . The low-latitude atmosphere also 
exh ibits scattered radiation of solar X -rays with a typical luminosity of ~ 3 x 10 15 erg 



s~ 1 (jBranduardi-Raymont et al. 2007bl ). Besides Jupiter itself, faint soft (< 1 keV) X-rays 
from the Io Plasma Torus (IPT) have been detected (~ 1 x 10 14 erg s _1 , Eisner et al. 2002). 
It may arise from bremsstrahlung emission by warm electrons with energies of few hundred 
to few thousand eV. Although we can expect hard X-rays (> 1 keV) from energetic keV to 
MeV electrons in the Jupiter's radiation belts, there has been no report on such an emission 
because past instruments lacked a necessary sensitivity to detect the extended hard X-ray 
emission. This indication motivated us to perform a decisive analysis using a deep 160 ks 
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observation data of Jupiter with Suzaku. 



Observations 



Suzaku (jMitsuda et al. 20071 ) o bserved Jupiter on 2 006 February 24-28 with the 
X-ray Imaging Spectrometer (XIS) (IKoyama et al. 20071 ). The XIS consists of three 
front-illuminated (FI) CCDs (XIS0, 2 and 3) and one back-illuminated (BI) CCD (XIS1). 

Due to the low-earth orbit of Suzaku and the large effective area, the XIS-FI has the 
lowest particle background among all X-ray CCDs in currently available X-ray observatories. 
Figure 5 in Mitsuda et al. (2007) summarizes the background normalized by the effective 
area and the field of view, which is a good measure of sensitivity determined by background 
for spatially extended emission. The background of the XIS FI is the lowest among X-ray 
CCDs onboard Suzaku, Chandra and XMM-Newton. It is 2 and 5 times lower than the 
XMM-Newton MOS and pn, and Chandra ACIS-I in 2-5 keV, respectively. For this reason, 
we use only the XIS data in this paper. 

During the observations, the spacecraft was repointed four times to allow for the 
planet's motion (0.7~1.8 arcmin per day). The XIS was operated in the normal mode. The 
data reduction was performed on the version 2.0.6.3 screened data provided by the Suzaku 
processing facility, using the HEAsoft analysis package ver 6.3.1. The net exposure of each 
FI and BI chip was 159 ks. For spectral fits, we generated response matrices and auxiliary 
files with xisrmfgen and xissimarf gen. 



Although an optical loading from Jupiter (visual magnitude my 



-2 during 



Suzaku observations) often becomes a problem in X-ray observations (jElsner et al. 20021 ) 
we concluded that it is negligible in this case, because optical loading effects such as 
degradation of the energy resolution, changes in the energy scales and/or the detection 
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efficiency were not seen in XIS spectra. For confirmation, we es timated the 



power of the XIS compared to CCDs onboard XMM-Newton ( ILumb 20001 ). We found 



optical blocking 



that the Suzaku XIS has almost the same blocking power as the XMM-Newton MOS 
with a thick filter, with which th e observation of Jupiter was successfully conducted 



( iBranduardi-Raymont et al. 2004] ) 



3. Extended X-ray Emission 

3.1. Imaging analysis 

To search for any emission from Jupiter, we have to remove Suzaku's orbital motion 
and transform the data into Jupiter's comoving frame using an ephemeris obtained from 
the Jet Propulsion Laboratory (JPL). Before conducting this procedure, we checked X-ray 
images without these corrections. We created two-band mosaic images using XIS1 (BI) and 
XISO+2+3 (FI) data as shown in figure HJ We corrected the data for positional difference in 
exposure times by using the exposure map generator xisexpmapgen. An extended emission 
is detected in both the 0.2-1 keV and 1-5 keV bands along the path of Jupiter. 

Beside the extended emission, we found signatures of emission from point sources near 
Jupiter in the 1-5 keV band. Since we found no known point sources in the X-ray source 
catalogue compiled by NASA_], we decided to detect sources from the Suzaku image. We 
utilized the wavelet function program wavdetect in the CIAO package Q We created four 
1-5 keV images with different bin sizes (1, 2, 4, and 8) and ran the point search program by 
setting the significant threshold above lxl0~ 6 , which corresponds to one spurious source 
in a 1024x1024 pixel map (XIS field of view). We summarized source lists and removed 

1 http: / /heasarc. gsfc.nasa.gov/cgi-bin/Tools/high_energy_source/high_energy _source.pl 
2 http:/ /cxc. harvard.edu/ciao/ 
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possible spurious detections by visual inspection. 

The obtained source positions are shown in figured] (b) (green circles with red lines). 
Twenty sources are detected with the 1-5 keV X-ray flux of 1 ~ 3 x 10~ 14 erg s _1 cm -2 . 
Here we converted a XIS count rate into flux, assuming a power-law spectrum with an 
photon index of 1.5 and an average absorption column toward this fielcjf , which represents a 
typical spectrum of the background active galac tic nucleus. This num ber is consistent with 
the canonical cosmic X-ray background model (jGiacconi et al. 20011 ). which predicts 27 
background sources above the detection limit. Although the source list contains marginally 
detected point sources, we decided to remove X-ray photons from the 1-5 keV image 
analysis afterward, for safety. We excluded circular regions with a diameter of 2 arcmin 
centered at individual sources. The diameter is equal to the beam size or the half power 
diameter of the XIS Q. Since all the sources are faint, contamination from the sources is 
thus negligible. Because there are no significant point sources around Jupiter in the 0.2-1 
keV image (fig. [Ha), we conducted this exclusion process only for the 1-5 keV image. 

We then corrected the two-band images for Suzaku's orbital motion and Jupiter's 
ephemeris. In case of the 1-5 keV image, to take into account the excluded regions, we 
created the exposure map for each observation and excluded the point source regions. Then 
we transformed the exposure map considering the Jupiter's orbit and divided the orbital 
motion corrected image by the corrected map. In case of the 0.2-1 keV image, we made the 
same analysis except for excluding the point source regions. 

The results obtained are shown in figure |5J We detected a significant hard X-ray 
emission extended over ~ 16 x 8 Rj which is positionally associated with the Jupiter's 



3 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 
4 http:/ /www. astro jp/ suzaku /doc/ suzaku_td 
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radiation belts, and Io's orbit. Because Io's orbit coincides with the IPT, the emission 
is also spatially associated with the IPT. The X-ray luminosity of the extended emission 
arises from (3.6±0.4) xlO 15 erg s" 1 and (3.3±0.5) xlO 15 erg s" 1 in 0.2-1 and 1-5 keV, 
respectively. Here and below the uncertainties are 68 % confidence intervals. We assumed a 
distance to Jupiter of 5.0 AU on 24-28 Feb 2006. An estimated contamination of the point 
sources outside the excluded regions to the extended emission is ~20 %. Hence, most of the 
emission is due to the extended component. 

To check the distribution of the emission, we created projection profiles along the 
horizontal axis as shown in figure [3J The hard X-ray emission is significantly extended 
more than the point spread function (green line). We found that the profile can be roughly 
represented by a simulated model using xissimarf gen, which is a sum of a point source 
and a elliptical uniform emission with radii of 4 and 1.3 arcmin (12 and 4 R$). The emission 
is thus extended over a wide region (> 6 On the other hand, the soft emission is 
marginally consistent with the point source. There is an excess on the right side of the soft 
X-ray peak, which suggests emission from the IPT that has been reported by Eisner et al. 
(2002). 

The hard X-rays are hence unlikely to originate only from a single point source or 
Jupiter. Furthermore, an artificial effect on the extended morphology due to the analysis 
can be excluded because of the different morphology from the soft X-ray image. The 
emission thus seems to contain truly diffuse component related to Jovian magnetospheric 
processes. 

We note that the faintness of the emission should have hindered its detection in 
the past observations. Considering the angular resolution of XMM-Newton, its surface 
brightness will be l/30~l/40 of Jupiter's aurora, From figure 5 in Branduardi-Raymont 
et al. (2007a), we can roughly estimate the CCD background level as 1/10 of the auroral 
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emission. Therefore the detection of the diffuse emission shuold be quite difficult. If the 
emission is related to Jovian radiation belts, there can be also time variation (see §4). 

3.2. Spectral analysis 

To investigate characteristics of the emission, we analyzed its X-ray spectrum by 
extracting photons around the extended emission region. For simplicity, we excluded the 
point sources from the event files without orbital motion corrections and extracted events 
from a circular region with a radius of 3 arcmin centered at the path of Jupiter. The 
surrounding region with an outer radius of 6 arcmin after subtracting the point sources 
was utilized as a background. Consequently, this extended emission region includes both 
Jupiter and the extended emission, because the limited angular resolution hindered us to 
spatially separate these two components. 

Figure S] shows the spectrum obtained. It consists of a flat continuum extending up 
to 5 keV and signs of lines around 0.25 keV and 0.56 keV. Hence, we fitted the spectrum 
with a power-law plus two Gaussian model. The best-fit model represents data well with 
X 2 jv ~0.5. No absorbing column was required, as expected from the source's proximity. 
The obtained photon index was 1.4±0.2, while the central energies of the two Gaussian 
were 0.24±0.1 keV and 0.56±0.1 keV. Fluxes for the 0.24 keV and 0.56 keV lines were 
(3.9±0.6) xl0~ 5 and (2.0±0.5) xl0~ 5 photons cm" 2 s _1 , respectively. 

Characteristics of the two lines such as a central energy and a flux were similar to 
those of the Jupiter's auroral emission observed with XMM-Newton in Nov 2003 (see 
table 3 in Branduardi-Raymont et al. 2007a). The two lines at 0.24 keV and 0.56 keV 
are presumably a line complex from C 5+ and ionized Mg, Si, S and a K Q emission from 
6+ , respectively. Such line emissions are distinct features of the charge exchange reaction 
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( ICravens et al. 20031 ). Taken together with the fact that the 0.2-1 keV X-rays arise from a 
compact region near Jupiter, most of soft X-rays likely originate from Jupiter's aurora via 
the charge exchange. 

In contrast, the 1-5 keV X-ray emission was represented by a simple flat continuum 
with a photon index of 1.4. Since the peak of the X-rays coincides with Jupiter (fig. [2]b), 
some fractions of the emission can arise from Jupiter. 



The power-law spectrum was in fact seen in the 



of Jupiter's aurora ( IBranduardi-Raymont et al. 20041 ). For comparison, we plotted 



past XMM-Newton observations 



auroral emission models in three past XMM-Newton observations (tables 2 and 3 in 
Branduardi-Raymont et al. 2007a) in figure HI The 1-5 keV fluxes of these models are 
smaller than that of the Suzaku best-fit model by a factor of 2~3. This is consistent with 
an independent estimation from the projection profile that most of the extended emission 
comes from diffuse emission other than Jupiter. 



4. Discussion 

We have discovered extended hard X-ray emission around Jupiter. Although its 
spatial association with Jupiter strongly suggests a diffuse component related to Jovian 
magnetospheric processes, we first have to examine serendipitous background sources. 
Because we corrected the image for the Jupiter's orbital motion, any faint background 
source on the Jupiter's path, that is missed in our point source identification procedure, can 
be seen as extended emission in the orbital motion corrected image. Considering the almost 
symmetric projection profile (fig. [3]) centered at Jupiter, the spatial distribution of the 
background sources should be symmetric about Jupiter. The most plausible case is that the 
emission is composed of Jupiter itself and a single background point source near the center 
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of the Jupiter's orbit. To test this hypothesis, we excluded the central part of the Jupiter's 
orbit with a circle region in the same way as figure [TJ (b), and created the orbital motion 
corrected image and projection profile in 1-5 keV. We still observed a significant extended 
emission around Jupiter, although the total count rate decreased by ~15 %. Hence, a single 
background source cannot fully explain the observed emission. Since multiple point sources 
with symmetric spatial distribution against Jupiter are less probable, below we discuss the 
origin of the emission, considering that the emission is truly diffuse. 

The flat power-law continuum of the emission suggests a non-thermal mechanism. There 
are three candidates for the emission mechanism; synchrotron emission, bremsstrahlung 
emission, and inverse-Compton scattering. The synchrotron interpretation can be easily 
rejected in terms of the required electron energy. From the X-ray image, we can observe 
the hard X-rays at the Io's orbital plane (a distance from Jupiter, r =5.9 Rj) or even more 
distant area. At r ~6 Rj, the magnetic field of the Jovian magnetosphere is quite weak on 
the order of 0.01 G. With such a weak magnetic field, if we want to generate synchrotron 
emission in the X-ray wavelength range, TeV electrons are necessary. This energy is more 
than 10 4 tim es higher than the observed maximum electron energy in the radiation belts 



(~50 MeV) (IBolton et al. 20021 ). Even if TeV electrons exist, their Larmor radius (3 x 10 11 
cm for 1 TeV electron at r =6 i?j = 4 x 10 10 cm) will be too large, so that TeV electrons 
are not trapped in the magnetic field line and will escape from the Jupiter's radiation belts. 

The bremsstrahlung emission is possible in terms of the electron energy, because 
we need only keV el ectrons which have been observed in the past in-situ measurements 
( jBagenal et al. 20041 ) . Such electrons will interact in an ambient medium and produce not 
only continuum but also X-ray lines by collisional excitation and inner-shell ionization (see 
§3.1 in Tatischeff 2002). From in-situ measurements, the plasma at r ~6 i?j is known to be 
mainly composed of heavy ions such as S + and + (see fig. 23.2 in Bagenal et al. 2004). 
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We estimated an equivalent width of S K a line (2.3 keV). We used cross sections for K-shell 
ionization by electron impact from Santos et al. (2003) and the bremsstrahlung loss rate 
of electrons given by Skibo et al. (1996), to estimate the line and continuum intensities, 
respectively. The calculated equivalent width of the S K a line was as large as ~10 keV, 
which should be observed in the XIS spectrum with the good energy resolution (~130 eV 
at 6 keV, Koyama et al. 2007). As shown in figure HJ no such line exists around 2 keV. 
Hence, the bremsstrahlung interpretation can be rejected. 

The remaining possibility is inverse- Compton scattering. We consider the scattering of 
solar light by ultra-relativistic electrons in the radiation belts. In this process, the scattered 
photon energy will be 



8 keV (Ep h /1.4 eV) (E c /50 MeV) 2 



(1) 



where E ph is an energy o 



(IRvbiki fc Lightman 19791 ) . 



solar optical photons and E e is that of relativistic electrons 
Considering the ma ximum observed electron energy of 50 MeV 



in the inner radiation belts (IBolton et al. 20021 ). the inverse- Compton scattering is possible 



in terms of the electron energy. Although 50 MeV electrons have not been observationally 
well confirmed at outer regions (r > several Rj), there may exist such ultra-relativistic 
electrons in these regions as expected from the empirical model (e.g., Divine & Garret 
1983). Also, since this mechanism does not need the emission line, the observed spectral 
shape can be safely explained. 

Assuming the inverse Compton scattering, we can estimate the index of the 
electron number density s pectrum as ~ 2, from the photon index of the X-ray spectrum 



(IRybiki fc Lightman 19791 ) . This is consistent wi th the Divine-Garrett (DG) empirical 



charged particle model of Jovian magneto sphere 



another model based on the Galileo data fjGarret et al. 20031 ) (~ 3 at 8 i?j). 



Divine fc Garret 19831) (~ 3 at 6 Rj) and 



We then proceeded to estimate the necessary electron number density. For simplicity, 
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we assumed that the emission region is an oblate spheroid with radii of 8, 8, and 4 
and the monotonic electron energy of 50 MeV. We here considered anisotropic angular 



distribution of the inverse- Compton scattering (IBrunetti et al. 20011 ). During the Suzaku 
observation, the Sun- Jupiter-Earth angle was 10°. Because of the larger cross section of the 
back scattering, the inverse- Compton flux will increase by a factor of ~3. The estimated 
average electron number density to explain 80% of the hard X-rays was ~0.005 cm -3 . 

For comparison, we utilized the DG empirical model, which provides a theoretical 
electron spectrum as a function of the distance from Jupiter. The DG model is based 
on the data taken with Pioneer and Voyager. To date, it is the best model to estimate 
the MeV electron distribution in the inner radiation belts (< 8 R s ). We found that 
the >20 MeV electron number density is only ~0.0007 and 0.0001 cm -3 at 4 and 6 Rj, 
respectively. The density will decrease as the distance increases. Hence, there is a factor of 
7~50 discrepan cy. Although the s ituation can be somewhat relaxed by an uncertainty of 



( IMivoshi et al. 19991 : 



the DG model (|Bolton et al. 200l[) and / or time variations of the trapped MeV electrons 



Bolton et al. 2002c ISantos-Costa et al. 20081 ). this interpretation 



needs further investigations. 

In conclusion, further X-ray, radio and in-situ observations are necessary to understand 
the origin of this diffuse hard X-ray emission, to identify the emission region with the 
Jupiter's magnetosphere, the IPT and/or else, and to clearly distinguish the emission from 
background point sources. 
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Fig. 1. — Suzaku XIS mosaic images of the vicinity of Jupiter in the (a) 0.2-1 keV (BI) 
and (b) 1-5 keV (FI) bands, displayed on the J2000.0 coordinates. Exposures are corrected 
and the count unit is counts s -1 binned pixel -1 . For clarity, images are binned by a factor 
of 8 and smoothed by a Gaussian of a = 5 pixels. Black circles show the size and posi- 
tion of Jupiter during the observations. Considering known pointing uncertainty of Suzaku 



(jUchiyama et al. 20081 ). positions of the circles are slightly shifted by +25 and —25 arcsec in 
the Ra and Dec directions, respectively, in order that the extended emission coincides with 
the Jupiter's path. Green circles with red lines mark excluded point source regions when we 
create the 1-5 keV image in comoving frame (see text). 
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Fig. 2. — Suzaku XIS images after correcting for the satellite's orbital motion and Jupiter's 
ephemeris in the (a) 0.2-1 keV (BI) and (b) 1-5 keV (FI) bands. The images are binned 
and smoothed in the same way as figured! In the panel (a), a circle indicates the expected 
position and size of Jupiter whose diameter is 39 arcsec. In the panel (b), grey lines indicate 
the equatorial crossing of magnetic field lines at 2, 4, 6, and 8 Jovian radius (Rj)- A black 
line is the path traced by Io. A photograph of Jupiter by Cassini taken from the JPL web 
site is overlaid. An arrow indicates the direction of the Sun. Boxes are used to obtain the 
projection profiles in figure [3j 
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Fig. 3. — Projection profiles along the horizontal axis, extracted from box regions in figure EJ 
Black and red points show the 0.2-1 and 1-5 keV data, respectively. Errors are la statistical 
ones. Green lines are a simplified model of Jupiter's X-rays (42 arcsec diameter) in the two 
energy bands, whose normalization and offset are tuned to coi ncide with the profile. The 



point spread function is almost the same at different energies (jSerlemitsos et al. 20071 ) . A 
dashed blue line indicates the expected profile of an uniform emission extended over an 
elliptical region with radii of 4 and 1.3 arcmin. A solid black line is the sum of the dashed 
green and blue lines. Dashed black lines are background levels estimated by fitting the data 
above the projected distance >14 arcmin with a constant. 
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Fig. 4. — Background subtracted BI (red) and FI (black) spectrum of the extended emission 
region, compared with the best-fit power-law plus two Gaussian models (solid line). Two 
Gaussians are shown in dashed and dash-dotted lines. Purple, green and blue dashed lines 
plotted for the FI spectrum are Jupiter's auroral continuum emission models in XMM-Newton 
observations (Branduardi-Raymont et al. 2007a). 



